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Processing signals from a set of photodetectors in gamma cameras - using cellular architecture 
so only those photodetectors involved in a given interaction 
are taken into account leading to high count rate 



The method for processing the signals produced by a set of photodetectors (60), in response to an 
event to be detected, uses 'a cellular architecture (88, 90-1... 90-6, 94) associated with each 
photodetector (60). The flux each photodetector receives can be compared with the flux of the 
neighbouring photodetectors (160-1. ..,160-6) to determine which set of photodetectors should be 
taken into account for a given interaction. Characteristic quantities, e.g., energy and position of the 
interaction, can then be calculated. 

- USE - In nuclear medicine e.g. in gamma cameras, positron emission tomography and correction of 
attenuation by transmission. 

- ADVANTAGE - High count rate, as only the signals from those photodetectors involved in a 
particular event are retained. Account can be taken of 2 events produced simultaneously at different 
places. (Dwg.7/15) 
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(54)pROCEDE ET DISPOSITIF DE TRAITEMENT DES SIGNAUX D'UN ENSEMBLE DE PHOTODETECTEURS 
W AYANT UNE ARCHITECTURE CELLULAIRE, ET APPLICATION AUX GAMMA-CAMERAS. 



(57) L'invention concerne un procede de traitement de si- 
grTaux produits par un ensemble de N photodetecteurs, en 
reponse a un evenement a detecter, dans lequel une archi- 
tecture cellulaire (88, 90-1, ( , 90-6, 92, 94), associee a 
chaque photodetecteur (60), permet de comparer le flux 
qu'il recoit avec les flux des photodetecteurs voisins (160- 
1, ( , 160-6), et de determiner quel est Pensemble de photo- 
detecteurs a considerer pour une interaction donnee. Des 
grandeurs caracteristiques de Interaction peuvent ensuite 
etre calculees. 
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PROCEDE ET DISPOSITIF DE TRAITEMENT DES SIGNAUX D 1 UN 
ENSEMBLE DE PHOTODETECTEURS AYANT UNE ARCHITECTURE 
CEliLUIiAIRE , ET APPLICATION AUX GAMMA- CAMERAS 

5 DESCRIPTION 

Domaine technique 

La presente invention concerne le traitement 
des donnees ou des signaux issus de, ou emis par, un 
ensemble de N photodetecteurs, en reponse a un 
10 evenement a caracteriser ou a identifier, par exemple a 
localiser . 

L f invention s 1 applique en particulier a la 
determination des caracteristiques, par exemple 
d'energie et de position, d'un evenement a partir de 

15 signaux fournis par des photomultiplicateurs equipant 
une gamma-camera ; par exemple, la position d r un 
evenement est reperee par rapport aux 

photomultiplicateurs eux-memes. On entend par gamma- 
camera une camera sensible au rayonnement gamma ( y ) - 

20 De telles cameras sont utilisees notamment a des fins 
d'imagerie medicale. 

Etat de la technique anterieure 

A I'heure actuelle, la plupart des gamma- 
25 cameras utilisees en medecine nucleaire sont des 

cameras fonctionnant selon le principe des cameras de 

type Anger. On peut se reporter a ce sujet au document 

US-3 Oil 057. 

Les gamma-cameras permettent en particulier de 
30 visualiser la repartition, dans un organe, de molecules 

marquees par • un isotope radioactif prealablement 

injecte au patient. 
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La structure et le f onct ionnement d ' une gamma- 
camera connue sont exposes et resumes ci-apres en 
reference aux figures 1, 2A et 2B annexees. 

La figure 1 montre une tete de detection 10 
5 d'une gamma-camera disposee en face d'un organe 12 
contenant des molecules marquees par un isotope radio- 
actif. 

La tete de detection 10 comporte un collimateur 
20, un cristal scintillateur 22, un guide de lumiere 24 
10 et une pluralite de tubes photomultiplicateurs 26 
juxtaposes de fagon a couvrir une face du guide de 
lumiere 24. Le scintillateur est, par exemple, un 

cristal de Nal (T £ ) . 

Le collimateur 20 a pour f onct ion de 

15 selectionner parmi tous les rayonnements gamma 30 emis 
par 1' organe 12 ceux qui atteignent la tete de 
detection sensiblement sous incidence normale. Le 
caractere selectif du collimateur permet d'augmenter la 
resolution et la nettete de 1 ' image produite. 

20 Cependant, 1 1 augmentation de la resolution se fait au 
detriment de la sensibilite. A titre d 1 exemple, pour 
environ 10000 photons y emis par 1' organe 12, un seul 
photon est ef f ectivement detecte. 

Les photons y ayant traverse le collimateur 

25 atteignent le cristal scintillateur 22 ou quasiment 
chaque " photon y est convert! en une pluralite de 
photons lumineux. Dans la suite du texte on designe par 
evenement chaque interaction d'un photon gamma avec le 
cristal, provoquant une scintillation. 

30 Les photomultiplicateurs 26 sont congus pour 

emettre une impulsion , electrique proportionnelle au 
nombre de photons lumineux regus du scintillateur pour 
chaque evenement . 
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Pour qu'un evenement de scintillation puisse 
etre localise plus precisement , les 

photomultiplicateurs 26 ne sont pas directement accoles 
au cristal scintillateur 22 mais sont separes de ce 
dernier par le guide de lumiere 24 . 

Les photomultiplicateurs emettent un signal 
dont 1' amplitude est proportionnelle a la quantite 
totale de lumiere produite dans le scintillateur par un 
rayonnement gamma, c'est-a-dire proportionnelle a son 
energie. Toutefois, le signal individuel de chaque 
photomultiplicateur depend aussi de la distance qui le 
separe du point d' interaction 30 du rayonnement gamma 
avec la matiere du scintillateur. En effet, chaque 
photomultiplicateur delivre une impulsion de courant 
proportionnelle au flux lumineux qu'il a recu. Dans 
l'exemple de la figure 1, des petits graphiques A, B, C 
montrent que des photomultiplicateurs 26a, 26b et 26c 
situes a differentes distances d'un point d ' interaction 
30 delivrent des signaux avec des amplitudes 
differentes . 

La position du point d ' interaction 30 d'un 
photon gamma est calculee dans la gamma-camera a partir 
des signaux provenant de 1 ' ensemble des 
photomultiplicateurs en effectuant une ponderation 
barycentrique des contributions de chaque 

photomultiplicateur . 

Le principe de la ponderation barycentrique tel 
qu'il est mis en oeuvre dans les cameras de type Anger 
apparait plus clairement en se reportant aux figures 2A 

et 2B annexees. 

La figure 2A montre le cablage electrique d'une 
tete de detection 10 d'une gamma-camera, qui relie 
cette camera a une unite de formation d'une image. La 
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tete de detection comporte une pluralite de 
photomultiplicateurs 26. 

Comme le montre la figure 2B, chaque 
photomultiplicateur 26 de la tete de detection est 
associe a quatre resistances notees RX~, RX + , RY~ et 
RY + . Les valeurs de ces resistances sont propres a 
chaque photomultiplicateur et dependent de la position 
du photomultiplicateur dans la tete de detection 10. 

Les resistances RX~, RX + , RY" et RY + de chaque 
photomultiplicateur sont reliees a la sortie 50 dudit 
photomultiplicateur, representee sur la figure 2B avec 
un symbole de generateur de courant. Elles sont d' autre 
part respectivement reliees a des lignes collectrices 
communes notees LX~, LX + , LY~, LY + , sur la figure 2A. 

Les lignes LX~, LX + , LY" et LY + sont a leur tour 
reliees respectivement a des integrateurs analogiques 
52X~, 52X + , 52Y~, 52Y + , et, par 1 ' intermediate de ceux- 
ci a des convertisseurs analogiques/numeriques 54X~, 
54X + , 54Y~, 54Y + . La sortie des convertisseurs 54X", 
54X + , 54 Y~, 54Y* est dirigee vers . un operateur numerique 
56. Les lignes LX~, LX + , LY~, LY + sont par ailleurs 
reliees a une voie commune, dite voie energie. Cette 
voie comporte egalement un integrateur 57 et un 
convertisseur analogique/numerique 58 et sa sortie est 
aussi dirigee vers 1 1 operateur 56. 

Grace au dispositif de la figure 2, on calcule 
\a position de 1 r interaction selon les equations 
suivantes (US-4 672 462) : 
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X+-X- 

x++x- 



et 



v=! Y+-Y- 
Y++Y" 



dans lesquelles X et Y indiquent les coordonnees selon 
deux directions orthogonales de la position de 
5 1 1 interaction sur le cristal et dans lesquelles X + , X", 
Y + , Y~ indiquent respectivement les signaux ponderes 
delivres par les integrateurs 52X + , 52X~, 52Y + , 52Y~. 

Les valeurs de X et Y de meme que I'energie 
totale E du rayon gamma ayant interagi , avec le cristal 

10 sont etablies par l'operateur numerique 56. Ces valeurs 
sont utilisees ensuite pour la construction d'une image 
comme decrit, par exemple, dans le document 
FR-2 669 439. 

Le calcul de la position de 1 1 interaction est 

15. entache d'une incertitude liee aux fluctuations 
statistiques de Poisson du nombre de photons lumineux 
et du nombre de photoelectrons produits pour chaque 
evenement, c'est-a-dire pour chaque photon gamma 
detecte. L'ecart type de la fluctuation est d'autant 

20 plus faible que le nombre de photons ou de 
photoelectrons est eleve. En raison de ce phenomene, il 
convient de collecter la lumiere le plus soigneusement 
possible. La resolution spatiale intrinseque de la 
camera est caracterisee par la largeur a mi-hauteur de 

25 la distribution des positions calculees pour une meme 
source ponctuelle collimatee posee sur le cristal 
scintillateur . 
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Pour des rayons gamma d'une energie de 140 keV, 
la resolution est generalement de 1'ordre de 3 a 4 mm. 

L' energie d'un photon gamma detecte est 
calculee en faisant la somme des contributions de tous 
les photomultiplicateurs ayant regu de la lumiere. Elle 
est aussi entachee d'une fluctuation statistique. La 
resolution en energie de la camera est caracterisee par 
le rapport de la largeur a mi-hauteur de la 
distribution des energies calculees a la valeur moyenne 
de la distribution, pour une meme source. 

La . resolution en energie est generalement de 
1'ordre de 9 a 11% pour des rayons gamma d'une energie 
de 14 0 keV. 

Finalement, une gamma-camera de type Anger 
presente I'avantage de permettre de calculer en temps 
reel le barycentre des signaux des photomultiplicateurs 
avec des moyens tres simples. 

En effet, le systeme decrit precedemment 
comporte un nombre limite de composants. De plus, les 
resistances utilisees pour injecter le signal des 
photomultiplicateurs dans les lignes collectrices sont 
tres peu couteuses. 

Une telle camera presente cependant aussi un 
desavantage majeur qui est un taux de comptage reduit. 
On entend par taux de comptage le nombre d'evenements, 
c'est-a-dire d 1 interactions entre un photon y et le 
scintillateur , que la camera est capable de traiter par 
unite de temps . 

Une des limitations du taux de comptage 
provient notamment du fait que la camera est incapable 
de traiter deux evenements ayant lieu sensiblement 
simultanement en des points distincts du cristal 
.scintillateur . 
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En effet, des evenements simultanes mais 
geometriquement distincts donnent naissance a des 
signaux electriques qui s'empilent dans les lignes 
collectrices LX~, LX + , LY" et LY + et qui ne peuvent plus 
etre distingues. Ces evenements sont aussi "perdus" 
pour la formation d 1 une image. 

Les gamma-cameras sont utilisees non seulement 
dans des. techniques d'imagerie traditionnelles, mais 
egalement dans deux autres techniques d'imagerie 
medicale ou la limitation du taux de comptage est aussi 
une contrainte redhibitoire . 

Ces techniques sont les techniques dites de 
"correction d ' attenuation par transmission" et de "PET 
(Positon Emission Tomography) en coincidence". 

La technique de correction d 1 attenuation par 
transmission consiste a tenir compte, lors de la 
formation d'une image medicale,, de 1 ' attenuation propre 
du tissu du patient entourant 1 1 organe examine. Pour 
connaitre cette attenuation, on mesure la transmission 
des rayonnements gamma vers une gamma-camera a travers 
le corps du patient. A cet effet on fait prendre place 
au patient entre une source externe tres active et la 
tete de detection de la gamma-camera. Ainsi, lors de la 
mesure du rayonnement transmis, un nombre eleve 
d' evenements ont lieu dans le cristal scintillateur . Le 
nombre . eleve d 1 evenements par unite de temps accroit 
aussi la probability d 1 avoir plusieurs evenements 
sensiblement simultanes. Une camera de type Anger 
classique se revele alors inappropriee . 

La technique de PET consiste a injecter au 
patient un element tel que F 18 apte a emettre des 
positons. L 1 annihilation d'un positon et d 1 un electron 
libere deux photons y emis dans des directions opposees 
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et ayant une energie de 511 keV. Ce phenomene physique 
est mis a profit dans la technique d'imagerie PET. Dans 
• cette technique on utilise une gamma-camera avec au 
moins deux tetes de detection disposees de part et 
d T autre du patient. Les tetes de detection utilisees ne 
sont pas equipees de collimateur. En effet, un 
traitement electronique des informations, dit 
traitement de coincidence, permet de selectionner parmi 
les evenements ceux qui coincident temporellement , et 
de calculer ainsi la trajectoire des photons gamma. 

Les tetes de detection sont done soumises a des 
flux eleves de rayonnement gamma. Les gamma-cameras 
classiques de type Anger ont un taux de comptage 
generalement trop limite pour une telle application. 

A titre indicatif, une gamma-camera de type 
Anger peut fonctionner normalement avec une detection 
de 1.10 5 evenements par seconde, tandis qu'en imagerie 
PET il faut au moins 1.10 6 evenements par seconde pour 
un fonctionnement normal. 

Une autre limitation des gamma-cameras du type 
Anger, decrites ci-dessus, tient au fait que le calcul 
du barycentre d'un evenement est fixe def initivement 
par la construction de la tete de detection et 

notamment par le choix des resistances RX~, RX + , RY~, 
RY + pour chaque photomul tiplicateur . De raeme, le calcul 

photomultiplicateurs sur une voie commune (voie 
energie) . 

Un dispositif classique attribue done un 
evenement a un photodetecteur et un seul, et prend en 
compte tous les signaux de tous les photodetecteurs 
pour le calcul de 1' energie. Or, les amplif icateurs des 
photomultiplicateurs non concernes par un evenement 
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ajoutent du bruit sur les lignes collectrices LX + , LX , 
LY + , LY~. 

Expose de 1' invention 

L 1 invention a pour objet un procede et un 
dispositif dans lesquels seule 1 f information, ou les 
signaux, issue des Ni photodetecteurs concernes, parmi 
un ensemble de N photodetecteurs, par un evenement, est 
numerisee et s.ubit un traitement local permettant d f en 
isoler le signal d'intensite maximum, ou d'amplitude 
maximale, et un certain nombre de signaux emis par les 
photodetecteurs voisins du photodetecteur ayant emis le 
signal d'intensite maximum, ou d'amplitude maximale. 

Plus precisement, l 1 invention a pour objet un 
procede de traitement de signaux produits par un 
ensemble de N photodetecteurs en reponse a un evenement 
a detecter, comportant : 

- la comparaison, avec une valeur seuil, pour chaque 
photodetecteur, d f une integrale numerique du signal 
numerise obtenu a partir du signal emis, 

- la comparaison, pour chaque photodetecteur emettant 
un signal dont I 1 integrale numerique depasse ledit 
seuil, de 1' integrale du signal qu'il emet avec 
1 ' integrale du signal emis avec chacun des 
photodetecteurs voisins, 

- la memorisation de la valeur de 1' integrale maximum 
parmi toutes les integrales comparees, ainsi que des 
integrales des signaux emis par les photodetecteurs 
voisins du photodetecteur dont le signal presente 
1' integrale maximum. 

II est ensuite possible de produire, a 1 1 aide 
des intensites memorisees, au moins un signal 
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representatif d'une grandeur caracteristique d'un 
evenement produit devant les photodetecteurs. 

Ainsi, seuls les signaux des photodetecteurs 
concernes par 1 1 interaction sont memorises ou 
5 contribuent au calcul, ou a la production d'un signal 
representatif, d'une grandeur caracteristique de 
1 1 evenement . II n'est pas necessaire de memoriser ou de 
tenir compte des signaux produits par tous les 
photodetecteurs, car le temps de lecture entrainerait 

10 • une forte degradation de la performance de taux de 
comptage. Chaque evenement ne concernant en principe 
qu'un petit nombre de photodetecteurs, le procede 
permet de ne retenir que ces photodetecteurs qui sont 
concernes par 1' evenement. De plus, deux evenements se 

15 produisant simultanement , ou quasi-simultanement , en 
deux endroits differents, par rapport a 1' ensemble de N 
photodetecteurs pourront alors etre pris en compte. 

L 1 invention a egalement pour objet un procede 
pour identifier une caracteristique d'un evenement se 

20 produisant devant un ensemble de N photodetecteurs, 
comportant : 

- 1 1 identification d'un sous-ensemble de Ni 
photodetecteurs, concernes par 1 ' interaction, c'est- 
a-dire ayant produit un signal en reponse a 

25 l f evenement, 

— la -production d'un signal representatif de la 
caracteristique de I 1 evenement, en fonction : 

* du signal d'intensite maximum parmi les N x signaux 
produits par les Ni photodetecteurs, en reponse a 

30 1' evenement, 

* et des signaux emis par les photodetecteurs 
voisins du photodetecteur dont le signal presente 
l'intensite maximum parmi les Ni photodetecteurs. 
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Ce second precede permet de generer des signaux 
representatif s d'une caracteristique d 1 un evenement se 
produisant devant 1' ensemble de N photodetecteurs , 
uniquement en fonction des signaux emis par les 
photodetecteurs concernes par ledit evenement. La 
encore, on evite une lecture inutile de l 1 ensemble du 
champ de N photodetecteurs, lecture qui consommerait 
beaucoup de temps et qui reduirait le taux de.comptage. 

L 1 invention a egalement pour objet des 
dispositifs permettant de mettre en oeuvre les procedes 
ci-dessus . 

L' invention a ete decrite en utilisant 
l 1 integrale numerique comme signal representatif de 
l'intensite de chaque signal provenant de chaque 
photodetecteur . 

En fait, on peut aussi utiliser, a la place de 
1 ' integrale numerique, l'intensite au maximum . du signal 
numerique ou du signal analogique. Cependant, 
1' integrale numerique est 1 1 information qui convient le 
mieux a la mise en oeuvre de I 1 invention. 



Breve description des figures 

De toute fagon, les caracteristiques et 
avantages de 1' invention apparaitront mieux a la 
lumiere de la description qui va suivre. Cette 
description porte sur les exemples de realisation, 
donnes a titre explicatif et non limitatif, en se 
referant a des dessins annexes sur lesquels : 

- la figure 1, deja decrite, est une coupe 
schematique d' une tete de detection d'une gamma-camera, 
par exemple de type Anger ; 

- les figures 2A et 2B, deja decrites, montrent 
de fagon schematique, un dispositif de l ! art anterieur, 
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pour la collection et le codage de signaux provenant de 
photomultiplicateurs de la tete de detection de la 
figure 1 ; 

- les figures 3A et 3B representent deux 
5 exemples de Champs de photodetecteurs , avec des 
symetries differentes ; 

\ - les figures 4, 5A et 5B representent une 
part\ie d'un dispositif associe a un photodetecteur ; 

\ - les figures 6A et 6B representent un signal 
10 analogique fourni par un photodetecteur (figure 6A) , 
ainsi \que le signal numerique correspondant (figure 
6B) ; \ 

\_ ia figure 7 represente un mode de realisation 

de .1 ' invention ; 
15 I ia figure 8 represente schematiquement un 

repere ulilise pour le calcul des coordonnees locales 
d'un evjenement par rapport a un champ de 
photodetecteurs ; 

- les figures 9A et 9B representent des signaux 
20 de photodetecteurs dans un cas d'evenement simple et 

dans un cas d'empilement ; 

-I la figure 10 est un chronogramme pour des 
signaux emis dans un dispositif selon 1' invention, 

- ^la figure 11 est un autre mode de realisation 

25 de 1' invention ; 

: - les figures 12 et 13 representent des 
'realisations d' elements fonctionnels dans un dispositif 

/ selon 1' invention ; 

- les figures 14 et 15 schematisent des 
'30 realisations 1 , de composants pour un dispositif selon 

w 1 ' invention . ' 
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Descri Etion djt ail lte de modes de realisation de 
1' invention 

L' exemple sera pris d'un ensemble de N 
photodetecteurs, ou de photomultiplicateurs , d'une 
gamma-camera. . Le nombre N est alors sensiblement 

compris entre 50 et 100. / 

Un tel ensemble de N photodetecteurs/ est de 
preference un reseau bi-dimensionnel de 

photodetecteurs, par exemple un reseau bi-dimensionnel 
de photomultiplicateurs. I 

Des exemples de reseaux bi-dimensionnels de 
photodetecteurs sent illustres sur les figures 3A et 
3B. Plus precisement, ces figures representent , dans 
chaque cas, une position d'un reseau de photodetecteurs 
en vue de dessus (ou, en d'autres termes, si.l'on se 
reporte a la figure 1, vue a partir du scintillateur 
22). La section transversale d'un photodetecteur peut 
avoir diverses formes, par exemple car/ree, hexagonale 

ou circulaire. 

La figure 3A represente, en vue de dessus, un 
champ mixte de photodetecteurs hexagonaAx et ronds . Les 
photodetecteurs a section hexagonale sbnt disposes en 

nid d'abeilles. J 

La figure 3B correspond a' un champ de 
photodetecteurs, en vue de flessus, chaque 
photodetecteur ayant une section carree. 

Lorsqu'un evenement se produi^ face au champ de 
photodetecteurs, cet evenement emettant un rayonnement 
dans une gamme spectrale correspondant a celle des 
photodetecteurs, un certain nombre / de photodetecteurs 
vont emettre des signaux. Sur la figure 3A, on a 
represente les ensembles de photodetecteurs concernes, 
e'est-a-dire emettant des signaux; pour deux evenements 
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se produisant a des endroits differentes par rapport au 
champ de photodetecteurs . Un premier evenement se 
produit face au photodetecteur 60 et celui-ci produit 
un signal. En reponse a cet evenement, les 
5 photodetecteurs voisins du photodetecteur 60 produisent 
egalement un signal : il s'agit, sur la figure 3A, des 
photodetecteurs 160 de premiere couronne (en traits 
hachures moyennement serres), immediatement voisins du 
photodetecteur 60, et des photodetecteurs 260 (en 

10 traits hachures espaces) , dits de deuxieme couronne 
autour du photodetecteur 60. Au-dela, les signaux emis 
par les autres photodetecteurs sont en general 
negligeables . De meme, un evenement se produisant face 
au photodetecteur 61 entraine une emission de signaux 

15 par le photodetecteur 61, mais aussi par les 
photodetecteurs de sa premiere couronne 161 (en traits 
hachures moyennement serres) et par les photodetecteurs 
261 de sa deuxieme couronne (en traits hachures 
espaces) . 

20 Sur la figure 3B, un evenement se produisant 

face a un photodetecteur 63 produit egalement un signal 
dans les photodetecteurs voisins 163 (de premiere 
couronne) et dans les photodetecteurs 263 (de deuxieme 
couronne) autour du photodetecteur 63. 

25 La figure 4 est une vue partielle d'un 

dispositif de traitement du signal d'un photodetecteur. 
Elle ne represente qu'un seul canal d'un tel 
dispositif, c'est-a-dire la partie du dispositif 
- associee a un unique photodetecteur 60. Des canaux 

30 similaires sont associes aux autres photodetecteurs, en 
particulier aux photodetecteurs de premiere et deuxieme 
couronnes qui entourent le photodetecteur 60. Le 
photodetecteur 60 est relie a un convert isseur courant- 
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tension 62. En reponse a un evenement detecte, on 
obtient sur la sortie 64 du convertisseur 62 un signal, 
par exemple du type de celui qui est illustre sur la 
figure 6A. 

5 Le graphique de la figure 6A indique, en 

ordonnees, 1' amplitude du signal correspondant a 
1' impulsion et, en abscisses, le temps. 1/ amplitude du 
signal et le temps sont indiques en echelles 
arbitraires. t 0 designe 1' instant de depart de 

10 1' impulsion fournie par le photodetecteur , et t x 
1' instant ou 1' impulsion redevient quasiment nulle, 
apres etre passee par un maximum. A titre indicatif, la 
duree correspondant a 1' intervalle ti - t 0 est de 
l'ordre d'une microseconde, dans le cas d'un 

15 photomultiplicateur d'une gamma-camera couple a un 

cristal de Nal (T£) . 

Le signal analogique present sur la borne de 
sortie 64 est dirige vers un convertisseur 
analogique-numerique 66. Ce dernier echantillonne 
20 chaque impulsion du signal en un certain nombre 
d' echantillons n, comme illustre sur la figure 6B. Deux 
echantillons consecutifs sont separes par un pas, ou 
intervalle d'horloge p (l'horloge fonctionnant a 
1/p Hz) . 

25 A titre d' exemple, le convertisseur 

. echantillonne chaque impulsion du signal en 
n = 10 echantillons. Pour un signal de 1 microseconde, 
un echantillonnage est alors effectue toutes les 100 
nanosecondes . 

30 Le convertisseur analogique-numerique 66 est, 

de preference, un convertisseur rapide, de type 
« flash » pouvant fonctionner a une frequence de 
l'ordre de 10 a 20 megahertz. 
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Le signal numerique issu du convertisseur 
analogique-numerique 66 dirige vers un integrateur 
numerique 68 (compose d'un registre a decalage a N 
bascules et d'un sommateur) . Cet integrateur effectue 
5 une somme glissante des echantillons qui lui sont 
transmis par le convertisseur analogique-numerique 66. 
La somme glissante est effectuee sur un nombre donne 
d' echantillons . Ce nombre predetermine est egal, par 
exemple, a 10. 

10 Le signal representatif de la valeur du signal 

code (en sortie de 66) est dirige vers des moyens 70 
aptes a detecter la presence d'un evenement . 

La figure 5A represente une realisation de la 
detection des evenements en detectant le passage du 

15 signal code par un maximum et en le validant lorsque ce 
signal est au-dessus d'un seuil predetermine. La sortie 
du codeur analogique numerique 66 est connectee a une 
des entrees d'un cornparateur 702. A la seconde entree 
701 est appliquee une valeur numerique predeterminee 

20 VP. Ce cornparateur- emet en sortie un signal 703 egal a 
0 si la valeur codee est inferieure a VP et 1 si la 
valeur codee est superieure a VP. Le signal 703 est 
connecte a une des entrees d'un ET 704. 

La sortie codee est aussi connectee a une des 

25 entrees 705 d'un deuxieme cornparateur 708 et a l r entree 
du registre 706 qui retarde la valeur codee d'un pas 
d'horloge. La sortie 707 de ce registre est connectee a 
la deuxieme entree 707 du cornparateur 708. Ce 
cornparateur emet une impulsion 709 egale a 1 lorsqu'il 

30 voit pour la premiere fois 1' entree 707 (qui represente 
la valeur codee a 1' instant t-1) superieure a 1 ' entree 
705 (qui represente la valeur codee a 1 ' instant t). La 
sortie 709 est egale a 0 dans les autres cas. Cette 
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15 



sortie 709 est connectee a la deuxieme entree du 
circuit ET, 704. Le signal 710 de sortie du ET 
represente 1 1 instant de passage du signal code par un 
maximum quand le signal code est superieur a VP. Le 
5 signal 710 est connecte a l'entree d'un circuit de 
decalage 712 qui le retarde d'un certain nombre de pas 
d'horloge. Ce nombre est regie de telle fagon que 
l r impulsion a sa sortie soit synchrone de 1' instant ou 
la valeur de sortie du sommateur glissant 68 est egale 
10 a l'integrale du signal detecte. 

La figure 5B represente une realisation du 
sommateur glissant 68 qui realise 1 ' integration de 
1 ' evenement . Le signal code en sortie du convertisseur 
analogique-numerique est dirige vers un registre a 
decalage 682 dont le nombre d'etages represente la 
duree d 1 integration . La sortie du registre a decalage 
682 est connectee a 1' entree soustractrice d 1 un 
soustracteur 684 tandis que 1 T entree positive est 
connectee au signal code. La sortie du soustracteur est 
envoyee sur un accumulateur 686. La sortie de 
1 1 aceumulateur 686 represente l'integrale de 
1' evenement lorsque le circuit 712 delivre son 
impulsion . 

L 1 ensemble des moyens, associe a un 
photodetecteur, et decrit ci-dessus en liaison avec la 
figure* 4 est designe par la reference 88.- A chaque 
photodetecteur est associe un tel ensemble de moyens. 

Chaque photodetecteur compare ensuite la valeur 
de son integrale numerique avec celle des 
photodetecteurs voisins. Un dispositif pour realiser 
cette comparaison est represente sur la figure 7. Sur 
cette figure, les photodetecteurs voisins, dits de 
premiere couronne, sont representes schematiquement par 



20 



25 
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les references 160-1 160-6. A chacun de ces 

photodetecteurs est associe un ensemble de moyens 
equivalent aux moyens 88 associes au photodetecteur 60. 
L'integrale numerique correspondant au signal delivre 
par le photodetecteur 60 est comparee avec l'integrale 
numerique du signal delivre par chacun des 

photodetecteurs 160-1 160-6, et ceci a l'aide de 

moyens de coraparaison, ou comparateurs 90-1, - .., 90-6. 
Les comparateurs travaillent en permanence (a chaque 
pas de codage) et fournissent un resultat de 
comparaison qui prend la valeur 1 lorsque la valeur de 
l'integrale courante du photodetecteur 60 est plus 
grande que l'integrale courante du photodetecteur 
voisin auquel il est compare, et la valeur 0 dans le 
cas contraire. Les 6 resultats de comparaison sont 
envoyes vers un dispositif 92 (figure 7) compose d'un 
circuit ET a 6 entrees, qui delivre un signal egal a 1 
lorsque tous les resultats de comparaison sont egaux a 
1, c'est-a-dire lorsque la valeur de l'integrale 
courante du photodetecteur 60 est plus grande que celle 
de tous ses voisins. La sortie du circuit ET est une 
premiere condition qui permet de valider, ou non, le 
signal EVT (ou demande de stockage DMDS) fourni par le 
circuit 712 (figure 5A) , et qui resulte de la detection 
d'un maximum sur le signal code. 

- Une fois ces valeurs memorisees, des moyens, 
Vi ' peuvent etre integres a 1' ensemble 94, peuvent 
permettre de calculer toute grandeur caracterist ique de 
1'evenement. Ainsi, il est possible de calculer 
1'energie totale de 1'evenement : pour cela, la somme 



E= ZEi de toutes les integrales numeriques E L 
i = 0 
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correspondant au photodetecteur principal 60 et aux 
photodetecteurs voisins 160-1, . 160-6 est realisee. 

II est egalement possible de calculer, ou 
d'estimer, les coordonnees locales de 1 1 evenement , par 
rapport au champ de N photodetecteurs, par exemple par 
des formules de type : 

6 

S EiXid 

X = (i) 

6 

£Ei 
i = 0 



Y = 



6 

i=0 



Y-; t 



6 

i=0 



(2) 



ou : 

- d represente la distance separant deux 
photodetecteurs adjacents sur 1 1 axe des X (entre les 
centres de ces deux photodetecteurs), 

- 1 est la distance separant deux photodetecteurs 
adjacents sur 1 1 axe des Y, 

- et ou les coordonnees des photodetecteurs 160-1, . -., 
160-6 sont reperees par leur centre sous la forme 
(X ±/ Yi) avec Xi€={-l,0,l} et Yie{-l,0,l}, dans un 
repere orthogonal (<D,i,j) tel qu'illustre sur la 
figure 8, ou O represente le milieu du photodetecteur 
central 60, i le vecteur unitaire suivant 1 ' axe des 
X, et j le vecteur unitaire suivant 1 ' axe des Y. 

La grandeur caracteristique de 1' evenement, une 
fois calculee, peut etre ensuite transmise a 
1 ' exterieur . 
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Parallelement aux comparaisons effectuees a 
I'aide des comparateurs 90-1, 90-6, la valeur de 

l'integrale numerique du photodetecteur 60 est 
transmise a un ensemble de comparateurs, dont chacun 
5 est associe a un des photodetecteur s peripheriques 
160-1, . . . , 160-6, le photodetecteur 60 etant alors 
traite comme photodetecteur de premiere couronne pour 
les autres photodetecteurs. Ainsi, a chaque 
photodetecteur de 1' ensemble du champ de N 

10 photodetecteurs est associe un dispositif du type de 
celui qui est illustre ci-dessus sur la figure 7, y 
compris un dispositif de memorisation 94 et, 
eventuellement de calcul . 

Du fait des comparaisons effectuees dans les 

15 comparateurs associes a chaque photodetecteur, c'est-a- 
dire du fait des comparaisons effectuees localement, 
c'est bien le photodetecteur produisant x un signal, 
d'intensite ou d'energie, maximale qui est considere 
comme photodetecteur "central" pour l'evenement, ses 

20 photodetecteurs voisins etant alors consideres comme 
photodetecteurs de premiere couronne pour 1 ' evenement ; 
de meme, c'est bien alors l'integrale numerique de ce 
photodetecteur central qui est stockee ou memorisee en 
tant que signal "principal", tandis que les integrales 

25 numeriques correspondant aux photodetecteurs de 
premiere couronne sont stockees ou memorisees dans les 
autres registres, et correspondent a des valeurs 
secondaires, chacune etant inferieure a la valeur de 
l'integrale numerique associee au photodetecteur 

30 central ou principal. 

Les procedes et dispositifs decrits ci-dessus 
sont ceux d'une "architecture cellulaire", qui permet 
aux photodetecteurs de comparer entre eux les flux 
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lumineux qu ? ils regoivent, et de determiner lesquels, 
parmi eux, sont concernes par 1 ' interaction . Un tel 
procede permet d'ameliorer la performance de taux de 
comptage des photodetecteurs, grace au traitement local 
5 d'une information numerisee. 

En outre, . deux evenements se . produisant 
simultanement, en deux endroits suffisamment distincts 
l f un de 1' autre, peuvent etre pris en compte. La 
encore, ceci permet d'accroitre le taux de comptage. 

10 En ce qui concerne le calcul de la, ou des, 

caracteristique (s) de l'evenement, une fois que le 
photodetecteur central et les photodetecteurs voisins 
ont ete identifies, deux possibilites sont disponibles . 

Selon la ■ premiere, les calculs des 

15 caracteristiques (par exemple : energie, position) sont 
effectues au niveau local, c'est-a-dire, dans le cas de 
la figure 7, a l'interieur meme du dispositif 94 
associe au photodetecteur 60. Les resultats des calculs 
sont ensuite diriges, par un bus commun a tous les 

20 photodetecteurs vers des moyens de collecte et 
d' exploitation exterieurs, par exemple un 

microprocesseur d ! un micro-ordinateur . C'est le procede 
qui a ete decrit ci-dessus. 

Selon une deuxieme methode, les calculs sont 

25 effectues a l'exterieur de l r ensemble des architectures 
cellulaires des photodetecteurs. Dans ce cas, toutes 
les valeurs memorisees dans les moyens de memorisation 
94 sont transferees vers^ un organe de calcul externe, 
non represents sur les figures. 

30 En termes de taux de comptage, c'est-a-dire de 

nombres d' evenements trajLtes par seconde, la premiere 
possibilite est avantageuse, du fait que seules les 
grandeurs caracteristiques associees a I'evenement sont 
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£ transferer vers I'exterieur, au lieu de toutes les 
integrales numeriques des photodetecteurs concernes. En 
outre, il y a alors un systeme de calcul par 
photodetecteur. 

La seconde possibility off re 1'avantage de 
n'etre pas figee dans les grandeurs a calculer, et de 
pouvoir modifier, pour une grandeur donnee, la methode 
de calcul. Cette seconde possibility, qui peut etre 
realisee avec un systeme programmable (DSP, 
microprocesseurs, operateurs dedies, etc..) est plus 
souple, d'autant qu'il est alors possible de regagner 
en taux de comptage en multipliant les moyens de 
transfert et de calcul. 

Des aspects particuliers de 1' invention 
permettent de tenir compte des problemes eventuels 
poses par les evenements saturants et/ou par les 
empilements . 

Les evenements saturants sont des evenements de 
tres grande energie, par exemple des rayons cosmiqiies, 
qui peuvent traverser le detecteur en provoquant en 
sortie une impulsion saturante dont la duree peut etre 
tres longue par rapport a celle des signaux utiles. 
Afin d'eviter ces problemes, 1 f evenement detecte est 
invalide des qu'un echantillon du signal digitalise 
correspondant atteint un certain . niveau dynamique, par 
exemple le niveau maximum (255 si cas de codage sur 8 
*bit) . 

Un empilement correspond a la coincidence de 
deux evenements simultanes ou quasi-simultanes 
(separes, dans le cas d ' une gamma-camera, par exemple 
par moins de 1 ]_is) et proches 1 1 un de 1' autre, c'est-a- 
dire espaces au plus par deux couronnes de 
photodetecteurs . 



B 12543. 3/PM 



2755816 



. 23 

La figure 9B donne 1* allure du signal resultant 
d'un empilement : les formes d'ondes s 1 additionnent . 

Considerons un evenement simple d 1 amplitude 
maximale (figure 9A) . Soit t 3 1 ' instant ou le signal 
passe au-dessus d'un seuil S et t 4 1 1 instant ou le 
signal repasse sous le seuil. Posons AT=t 4 -t 3 . 

Considerons maintenant un evenement quelconque. 
Soit t 3 f 1' instant ou le signal passe au-dessus du 
seuil. L' evenement sera considere comme valide si, a 
I 1 instant t 3 f +AT, le signal est au seuil S ou en- 
dessous du seuil S. Ceci sera vrai pour un evenement 
simple d' amplitude moindre ; par contre en cas 
d 1 empilement (figure 9B) , le signal ne repassera au 
niveau S qu'a t 3 '+AT' et sera au-dessus du seuil a 
t 3 '+AT, puisque AT 1 >AT Dans un tel cas, egalement, 
1 ' evenement detecte est invalide. 

Pour eviter tous ces problemes de saturation ou 
d r empilement s, on repere 1 1 instant de f ranchissement du 
seuil S par le signal. Des le seuil franchi, un 
compteur se declenche. Si n coups d'horloge plus tard, . 
n etant fixe et correspondant a la duree AT, 
1 1 echantillon est superieur au seuil, alors 1 1 evenement 
est invalide : un signal de validite (VAL) prend la 
valeur r, 0 u . AT et S sont fixes en fonction du type 
d 1 evenement a detecter. 

On peut choisir que, si le signal d'un 
photodetecteur est invalide, il le reste pendant un 
certain temps, par exemple 1 ps. 

II a ete explique ci-dessus comment tenir 
compte, pour le calcul d'une caracteristique d'un 
evenement, des photodetecteurs de premiere couronne 
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(photodetecteurs 160, 161 sur la figure 3A) autour du 
photodetecteur principal. 

II peut etre important de prendre egalement en 
compte les photodetecteurs de deuxieme couronne . 

A cette fin, des qu 1 un photodetecteur est 
retenu en tant que photodetecteur central (c'est-a-dire 
que son signal est le plus grand et que, si un test de 
validite est realise, il est valide) , les moyens 92 
emettent un signal "wevt" vers le systeme de stockage 
local 94 et vers les systemes de stockage des 
photodetecteurs de premiere couronne indiquant qu'il 
faut retenir 1* evenement, c'est-a-dire qu'il faut 
stocker les valeurs des integrales et eventuellement 
calculer les grandeurs caracterist iques de 1 1 evenement . 
Un tel signal est envoye de maniere synchrone a un coup 
d 1 horloge . 

Comme, par ailleurs, un signal d 1 evenement 
"evt" est emis sur le coup d ! horloge precedent, le 
schema temporel resultant est celui de la figure 10. 

Apres que le signal "wevt" ait ete transmis aux 
voisins de la premiere couronne, ceux-ci transmettent 
la valeur de deux integrales de photodetecteurs de 
deuxieme couronne a la memoire 94 du photodetecteur 
central (figure 7) . 

Une fois memorisees les valeurs des integrales 
numeriques des photodetecteurs de premiere et de 
deuxieme couronnes, ces valeurs peuvent etre transmises 
vers l'exterieur, pour un calcul exterieur de grandeurs 
du type energie ou position de l f evenement. 

Ou bien, conformement a ce qui a deja ete 
explique ci-dessus, les grandeurs sont calculees par 
les moyens 94 eux-memes et ensuite transmises a 
1 1 exterieur . 
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L'une et 1 1 autre de ces possibilities presentent 
chacune des avantages qui ont deja ete indiques 
precedemment . 

Le calcul de 1 1 energie se fera par exemple en 
appliquant la formule : 

19 

E = Z E± 
i = l 

(ou 1=1, ...19, designe les 19 valeurs d 1 integrales 
numeriques a prendre en compte pour un systeme a 
structure hexagonale) . 

Le calcul de la position se fera par exemple en 
appliquant les formules suivantes, semblables aux 
formules (1) et (2) deja presentees : 

19 

X *iE±d 
X = ^ (3) 

ZEi 
i = l 

19 

z w 

Y - ±=± (4) 

ZEi 
i = l 

Xi et Y ± reperent les coordonnees des photodetecteurs de 
premiere et de deuxieme couronnes dans un repere du 
type de celui de la figure 8 : X i e{-2 / -l / 0,1, 2} ; 
¥^{-4,-3,-2,-1,0,1,2,3,4}. 

La figure 11 represente un autre mode de 
realisation d ' un dispositif electronique 170 associe a 
un photodetecteur 60. Un tel dispositif est associe par 
ailleurs a chacun des autres photodetecteurs. 
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Une horloge unique est appliquee a tous les 
photodetecteurs et definit le pas p de codage. 

Un tel dispositif comporte d'abord un premier 
ensemble fonctionnel 188 qui realise les fonctions 
suivantes, deja decrites ci-dessus : conversion 
courant-tension, digitalisation du signal issu du 
photodetecteur , integration, detection du passage du 
signal par un maximum. De plus, un signal DMDS 
(equivalent au signal EVT du mode de realisation 
precedent) de demande de stockage est genere a chaque 
detection d'un maximum et synchronise avec l f integrale 
du signal concerne disponible. 

En outre, cet ensemble fonctionnel peut assurer 
les fonctions de test de vraisemblance de forme de 
1' impulsion par analyse de 1' amplitude ou de la duree 
au-dessus d'un seuil, comme explique plus haut. 

Lors de chaque pas de codage, le dispositif 170 
transmet aux dispositifs correspondant s des 
photodetecteurs 160-1, 160-6 voisins les 'donnees 

suivantes : 

1 - la valeur de 1 1 integrale courante DATA 0 , 

2 - un signal de validite (signal VAL 0 ) emis lorsque les 

deux conditions suivantes sont remplies : 

* le dispositif 17 0 est apte a recevoir un ordre 
de stockage (registre de stockage 194 disponible) , 

* ■ le signal DATA 0 qu'il transmet a ses voisins 
est correct (signal teste VRAI) T 

3 - un signal de decision de stockage (signal DS 0 ) dont 

1 1 elaboration est expliquee par la suite. 

En retour, le dispositif 170 regoit done de 
chacun des voisins 106-1, . .., 160-6 de premiere 
couronne : 

1 -1' integrale courante DATAi (i=l,' 6), 
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2 - un signal de validite VALi (i=l, - . . , 6) , 

3 - un signal de decision de stockage DS± (i=l, 6), 

equivalent au signal WEVT du mode de realisation 
precedent . 

5 Ainsi, a chaque pas de codage, le dispositif. 

170 associe au photodetecteur 60 peut comparer la 
valeur courante DATA 0 avec celle de ses voisins et se 
determiner comme etant le plus grand ou non. Lorsqu'a 
un instant donne un eveneraent est detecte (passage du 

10 signal code par un maximum, superieur a un seuil) , et 
lorsque l'integrale numerique de cet evenement est 
disponible, le dispositif fait une demande de 
stockage : un signal DMDS est emis vers un dispositif 
de generation de decision de stockage 192 (equivalent 

15 au dispositif 92 deja decrit ci-dessus) . 

Cette demande est confirmee si, a cet instant : 

1 - la valeur de l'integrale du photodetecteur 60, qui 

demande le stockage, est plus grande que la valeur 
de l'integrale de chaque voisin de premiere 
20 couronne, 

2 - le critere de validite (signaux VAL 0 et VALi, i = l-6) 

est rempli pour le photodetecteur demandeur ainsi 
que pour ses voisins de premiere couronne. 

Lorsque ces conditions sont reunies, la demande 
25 de stockage est validee (signal DS 0 ) . Elle est alors 
effective dans le photodetecteur demandeur et transmise 
a ses voisins de premiere couronne par 1 ? intermediaire 
du signal "decision de stockage" DS 0 . 
Ceci a pour effet : 
30 1 - dans le dispositif 170, de stocker en memoire 194, 
la valeur DATA 0 de l'integrale courante du 
photodetecteur 60 et les integrales DATAi (i=l-6) de 
ses voisins de premiere couronne, 
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2 - dans chaque dispositif electronique associe a chacun 
des photodetecteurs voisins (de premiere couronne) 
de stocker la valeur courante des integrales de deux 
de ses propres photodetecteurs voisins de premiere 
couronne (qui sent, pour certains, des 
photodetecteurs de deuxieme couronne du 

photodetecteur demandeur) choisis de telle fagon que 
toutes les valeurs courantes des photodetecteurs de 
deuxieme couronne du photodetecteur demandeur soient 
recuperees . 

Lorsque cette operation est effectuee, le 
systeme se trouve dans I'etat suivant : 

- le registre de stockage 194 du photodetecteur 
demandeur 60 dispose de la valeur de sa propre 
integrale numerique et de la valeur de 1 1 integrale 
des photodetecteurs voisins de premiere couronne, 

- les registres des photodetecteurs voisins de premiere 
couronne disposent chacun de deux valeurs d ? integrale 
de photodetecteurs de deuxieme couronne. 

On peut preciser de la maniere suivante le role 
des entrees DSi sur le dispositif 194 : 

- sur le photodetecteur demandeur il ne peut pas y 
avoir de DSi actifs au moment ou le photodetecteur 
demandeur envoie une decision de stockage car les 
integrales des photodetecteurs de premiere couronne 
sont- plus petites que celles du photodetecteur 
central , 

- sur les photodetecteurs de premiere couronne, DSi 
sert a stocker la valeur de deux integrales de 
deuxieme couronne. 

Ainsi pour des photodetecteurs hexagonaux, ou 
ronds, arranges en nid d'abeilles (cas de la figure 
3A) , le registre de stockage du photodetecteur 
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demandeur 60 aura stocke sept valeurs ; chacun de ses 
six photodetecteurs de premiere couronne stocke deux 
valeurs, soit en tout dix-neuf valeurs. De meme pour 
des photodetecteurs carres (figure 3B) , le 
5 photodetecteur demandeur aura stocke neuf valeurs, et 
chacun de ses huit voisins deux valeurs, soit en tout 
vingt-cinq valeurs. 

Ceci est vrai pour les photodetecteurs 
demandeurs ayant une premiere et une deuxieme couronnes 

10 completes. Pour garder la modularite du circuit, 
lorsque les photodetecteurs n'ont pas tous leurs 
voisins (sur les bords du champ) , les entrees des 
photodetecteurs manquants seront fixees a zero. 

Des moyens 196 de type sequenceur permettent de 

15 realiser la gestion des transfert et du calcul et de 
donner aux moyens 198, de generation d'un signal 
d 1 occupation, un signal de fin d 1 utilisation du 
stockage . 

Les moyens de multiplexage 171 sont pilotes par 
20 la fonction 196 "gestion des transferts" et assurent la 
fonction suivante. Lorsqu'une decision de stockage a 
ete effectuee les photodetecteurs de premiere couronne 
ont memorise la valeur de deux photodetecteurs de 
deuxieme couronne. Pour transferer ces valeurs vers le 
25 photodetecteur demandeur on va utiliser, pour diminuer 
le nombre de liaisons entre photodetecteurs, le bus de 
sortie DATA de chaque photodetecteur de premiere 
couronne et remplacer la valeur de 1' integrale courante 
successivement par les deux valeurs a transmettre. Le 
30 photodetecteur demandeur va ainsi recevoir sur des 
entrees DATA a a DATA f les valeurs des photodetecteurs de 
deuxieme couronne. II les stocke a cote de celles qu'il 
avait deja. 
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Les moyens 192 de generation de decision de 
stockage peuvent etre realises de la maniere decrite en 
figure 12. 

Les . moyens 198 de generation d'un signal 
"occupation" peuvent avoir la forme illustree en figure 
13, ou la reference 193 designe une bascule de type D 
avec preset. 

La sortie occupation est raise a l'etat logique 
(VRAI) par -l 1 une des sept entrees decision de stockage. 
Cet etat est maintenu jusqu'a ce que les valeurs 
acquises aient ete utilises. Le systeme de gestion des 
transferts fournit alors une impulsion de remise a 
zero. 

Les fonctions et le systeme decrits ci-dessus 
peuvent etre integres dans un circuit ASIC. Autrement 
dit, a chaque photodetecteur du reseau de N 
photodetecteurs, est associe un circuit ASIC. 

L'ASIC numerique de chaque photodetecteur peut 
comporter ou assurer les fonctions suivantes : 
integration, detection du maximum, entrees des 
photodetecteurs voisins, sortie vers les 

photodetecteurs voisins, comparateurs , registres de 
stockage, sequencement , dispositif de calcul des ou 
d'une grandeur (s) caracteristique ( s ) d'un evenement, 
dispositif de lecture . 

L'ASIC numerique peut en outre comporter les 
fonctions suivantes : 

— multiplication de gain, 

— correction de ligne de base, 

— gestion des evenements saturants. 

Afin d'evaluer la cornplexite de l'ASIC 
numerique associe a chaque photodetecteur, c T est-a-dire 
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son nombre de portes elementaires, sa surface et son 

nombre de plots d ' entrees/sorties , on peut etudier ses 

composants qui requierent le plus de portes. Ce sont 

principalement : 

o les registres, 

o les additionneurs, 

o les compteurs, 

o les multiplexeurs , 

o les multiplieurs et diviseurs (pour le barycentre) . 

Tous ces composants sont constitues de cellules 
de base que sont : 
o la bascule D, 

» 1 'additionneur 1 bit (cellule FA detaillee plus 
loin) , 

o le multiplexeur 2 bits vers 1 bit. 

Ce a quoi il faut ajouter la logique 
combinatoire environnante . 

La bibliotheque MHS donne le nombre de portes 
logiques equivalent a chacune de ces cellules.de base. 
Ainsi : 

© une bascule D simple horloge ("Single-clock D flip 

flop" : DFF) comprend 5 portes, 
o un additionneur 1 bit (FA) comprend 9 porte, 
o un multiplexeur 2 vers 1 (MUX2) comprend 3 portes. 

Un additionneur de 2 nombres de n bits 
necessite n cellules FA ("Full Adder", constitute de 9 
portes ) . 

Du point de vue du nombre de portes, un 
soustracteur est a considerer comme equivalent a un 
additionneur. 

Pour multiplier un nombre de n bits par un 
nombre de m bits, avec n>m, un multiplieur comprend 
n*(m--l) cellules FA. A titre d'exemple, un multiplieur 
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parallele 13x13 bits de la bibliotheque ES2 a un temps 
de calcul de 3 6.8 ns pour une surface occupee de 
0,78 mm 2 , ce qui est amplement suffisant pour les 
performances envisagees ici. 
5 Le schema de la figure 14 donne un apercpu 

schematique de la fagon dont sont agencees les entites 
d 1 integration et de commande (les moyens 188 de la 
figure 12) . . 

Le sommateur glissant est detaille sur la 

10 figure 14, la partie assurant le multiplexage de la. 
sortie est detaillee sur la figure 15. 

Le sommateur comporte un registre a decalage a 
N bascules D x , D N et un sommateur. Le soustracteur 

220 permet de calculer l'ecart entre l f entree et la 

15 sortie du registre a decalage. 

L 1 additionneur 221 accumule ces ecarts dans le 
registre 223. Un circuit de detection de ligne de base 
222 compose d'un comparateur et d'un registre a 
decalage pilote un registre de stockage 225 De la ligne 

20 de base. Celle-ci est alors soustraite par le 
soustracteur 226 a la valeur courante de la somme 
glissante. Les moyens 224 qui assurent la detection des 
evenements et la notion de validite .sont 
essentiellement composes d'une detection de maximum, 

25 realisee a l'aide d'un comparateur, suivi d'un registre 
a decalage necessaire pour synchroniser le signal 
evenement avec la sortie de 1 1 integrale correspondante . 

Le sommateur glissant est compose d'un registre 
a decalage (10 registres a 8 bits soit 10*8=80 bascules 

30 D),d'un registre de 12 bits (12 bascules D) , d'un 
sommateur sur 12 bits et d'un soustracteur sur 8 bits, 
soit (12+8)=20 cellules FA. 

Le sommateur comprend done 92 DFF et 20 FA . 
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- La determination de la ligne de base requiert 
une logique de compa.raison et un registre a decalage de 
16 bits avec remise a zero, soit 16 DFF. 

II faut rajouter le registre d 1 offset 12 bits 
5 et le soustracteur 12 bits, ce qui correspond a 12 DFF 
et 12 FA. 

La partie multiplexage de la sortie peut se 
schematiser comme illustre en figure 15. 

Elle comporte 2 registres de 12 bits, (24 DFF), 
10 1 multiplexeur a 6 entrees de 12 bits, et 1 
multiplexeur a 8 entrees de 12 bits. 

Or, pour faire un MUX a 6 entrees, on a besoin 
de 5 MUX a 2 entrees. De meme, il faut 7 MUX a 2 
entrees pour faire un MUX a 8 entrees. 
15 Done au total il faut <7+8)*12=180 MUX2 . 

En resume, I'entite d 1 integration 1 et de 
multiplexage comprend : 
o 144 DFF, soit 144*5=720 portes, 
o 32 FA, soit 32*9=288 portes, 
20 o 180 MUX2, soit 180*3=540 portes. 

Ce qui porte le nombre de portes a un total 
d' environ 1600 portes, en comptant la logique 
combinatoire environnante . 

25 L'entite de commande (92, 192 sur les figures 7 

et 12) -comprend, notamment : 

4 comparateurs 8 bits (soit environ 4*8*5=160 
portes) . 

o Un compteur de 10 bits et un de 6 bits, soit 16 DFF, 
30 done 80 portes, 

«»* • si l 1 on rajoute toute la logique de sequencement , on 
peut evaluer l'entite commande a 300 portes. 
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L'entite de comparaison (190-1, 190-6 sur 

la figure 11) comporte : 

• 6 comparateurs 12 bits, soit 360 portes, 

• 19 registres de 12 bits pour le stockage des 
energies, soit 228 FFD dont 1200 portes. 

En ce qui concerne l'entite de calcul du 
barycentre, on a deja vu que ce calcul peut se faire en 
appliquant les formules (1), (2) et (3), (4) ci-dessus. 

La structure de calcul comprend done : 

• un multiplieur 12*4 bits, soit 36 FA, 

• un registre 14 bits, soit 14 DFF, 

• un sommateur 14 bits, soit 14 Fa, 

• il faut 3 autres sommateurs 14 bits pour obtenir la 
somme des 18 valeurs, soit 42 FA, 

• un registre de stockage du resultat sur 14 bits (X 1 
ou Y f ) , soit 14 DFF. 

Cette structure est a repeter deux fois 
puisqu'il y a deux coordonnees, done au total on a : 

• 384 FA, soit 3460 portes, 

• 112 DFF, soit 560 portes, 
done 4000 portes. 

On a ainsi besoin de deux multiplieurs 
paralleles, un 14*13 bits et un 14*6 bits, ce qui 
represente 238 FA. 

- II faut un registre de 27 bits et un de 20 bits 
pour stocker les resultats, soit 47 DFF. Cela donne au 
total 2400 portes. 

II reste a diviser ces nombres par I'energie 
totale, pour obtenir les coordonnees finales. Pour 
cela, on. peut se contenter d'effectuer une division 
sequentielle en 9 coups d'horloge. A cette fin, on peut 
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evaluer que, pour diviser 27 bits par 14 bits et 

obtenir 9 bits, on a besoin de : 

o 1 registre 27 bits avec 27 MUX2, 

o 1 registre 14 bits, 

o 1 registre 9 bits, 

o 1 soustracteur 14 bits. 

Ceci equivaut a 50 DFF, 14 FA et 27 MUX2 , soit 
au total environ 500 portes. 

II faut deux diviseurs au total, ce qui fait 
1000 portes. 

L'entite de "transf ert " 196 (figure 11) 
necessite trois registres pour stocker X, Y et E, soit 
deux registres de 9 bits et 1 de 14 bits, ce qui 
equivaut a 32 DFF, soit 160 portes. Si 1 1 on ajoute la 
logique de sequencement , on peut evaluer l'entite a 200 
portes . 

L'ASIC est ainsi evalue a 11000 portes. Dans la 
technologie de 0,7 yam, le nombre de portes par mm 2 est 
d' environ 1400 portes/mm 2 . On peut done evaluer la 
surface occupee par les portes a 8 mm2, soit environ 
3 mmx3 mm. 

Qn boitier de 140 pattes peut etre utilise. 

L 1 invention decrite s 1 applique a une gamma- 
camera telle que decrite dans 1 1 introduction a la 
demande, en liaison avec la figure 1. Elle s' applique 
avantageusement a la mise en oeuvre des techniques de 
"correction d 1 attenuation par transmission" et de "PET 
(Position Emission Tomography) en coincidence, deja 
decrites dans 1 1 introduction a la presente demande. 

L 1 invention a ete decrite en utilisant 
l'integrale numerique comme signal representatif de 
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l 1 intensite de chaque signal provenant de chaque 
photodetecteur . 

En fait, on peut aussi utiliser, a la place de 
1 ' integrale numerique, .1' intensite au maximum du signal 
numerique ou du signal analogique. Cependant, 
1 1 integrale numerique est la quantite qui convient le 
mieux a la mise en oeuvre de 1' invention. 

Si on veut mettre en oeuvre un dispositif dans 
lequel 1' intensite au maximum du signal numerique ou du 
signal analogique est . choisie eh tant que quantite 
representative de 1 ' intensite de chaque signal 
provenant de chaque photodetecteur , il suffit de 
remplacer, dans la chaine de detection associee a 
chaque photodetecteur, les moyens pour calculer une 
integrale numerique par des moyens pour detecter un 
maximum numerique ou analogique. La gestion des donnees 
obtenues (comparaison, memorisation) se fait ensuite de 
la meme maniere que decrite ci-dessus. 
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RE VEN PI CAT I ON 5 

1. Procede de traitement de signaux produits 
par un ensemble de N photodetecteurs en reponse a un 
evenement a detecter, comportant : 

- la comparaison, pour chaque photodetecteur, d T un 
signal representatif de I'intensite du signal 
numerise obtenu a partir du signal emis, avec une 
valeur seuil, 

- la comparaison, pour chaque photodetecteur emettant 
un signal dont I'intensite depasse ledit seuil, de 
I'intensite du signal qu'il emet avec I'intensite du 
signal emis par chacun des photodetecteurs voisins, 
dits de premiere couronne, 

- la memorisation de la valeur de I'intensite maximum 
parmi toutes les intensites comparees, ainsi que des 
intensites des signaux emis par les photodetecteurs 
voisins du photodetecteur dont le signal presente 
I'intensite maximum. 

2. Procede selon la revendication 1, le signal 
representatif de I'intensite du signal numerise etant 
soit une integrale numerique du signal numerise, soit 
I'intensite au maximum du signal numerise, soit 
I'intensite au maximum du signal analogique. 

3. Procede selon la revendication 1 ou 2, au 
moins un signal, representatif d'une grandeur 
caracteristique de 1' evenement, etant produit a 1 ' aide 
"cles intensites memorisees. 

4. Procede selon la revendication 3, un signal 
representatif de l'energie de 1' evenement etant 
produit, 

5. Procede selon la revendication 3, un signal 
representatif de la position de 1' evenement par rapport 
a 1' ensemble des N photodetecteurs etant produit. 
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6. Precede selon la revendicat ion 3, au moins 
deux signaux representatif des coordonnees 

barycentriques de 1 1 evenement par rapport a l 1 ensemble 

des N photodetecteurs etant produit. 

5 7. Precede ■ selon 1 1 une des revendications 1 a 

6, un signal d ' evenement etant emis, pour chaque 
photodetecteur , lorsque le signal representatif de 
l'intensite du signal numerique est forme. 

8. Precede selon la revendication 7, 
10 l'intensite du signal numerique etant representee par 

son integrale numerique, celle-ci etant calculee sur un 
certain nombre de pas d'horloge, a partir du maximum du 
signal numerique, et le signal d 1 evenement etant emis 
un certain nombre de pas d'horloge apres le passage par 
15 le maximum. . 

9. Precede selon 1 1 une des revendications 1 a 

8, comportant en outre une etape de detection et de 
rejet d'evenements saturants et/ou d 1 empilement s 
d 1 evenements . 

20 10. Procede selon 1 ' une des revendications 1 a 

9, comportant en outre une etape de calcul et de 
memorisation des intensites des photodetecteurs dits de 
deuxieme couronne . 

11. Procede selon 1 ' une des revendications 3 a 
25 6 et selon la revendication 10, le signal representatif 

d'une -grandeur caracteristique de 1 ' evenement etant 
produit a l'aide des intensites memorisees des 
photodetecteurs de premiere et deuxieme couronnes . 

12. Procede pour identifier une caracteristique 
30 d'un evenement se produisant devant un ensemble de N 

photodetecteurs, comportant : 

- 1 1 identification d'un sous-ensemble de Ni 
photodetecteurs, concernes par 1 1 interaction, e'est- 
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a-dire ayant produit un signal en reponse a 
l'evenement, 

- la production d 1 un signal representatif de la 
caracteristique 'de l'evenement, en fonction des 
intensites : 

* du signal d'intensite maximum parmi les Ni signaux 
produits par les Ni photodetecteurs , en reponse a 
l'evenement, 

* et des signaux emis par les photodetecteurs 
voisins (dits de premiere couronne) du 
photodetecteur dont le signal presente l'intensite 
maximum parmi les Ni photodetecteurs. 

13. Procede selon la revendication 12, la 
caracteristique etant I'energie de l'evenement. 

14. Procede selon la revendication 12, la 
caracteristique etant la position de l'evenement par 
rapport a 1' ensemble de N photodetecteurs. 

15. Procede selon la , revendication 14, la 
position etant reperee par un couple de coordonnees 
barycentriques . 

16. Procede selon la revendication 12, le 
signal representatif de la caracteristique d'un 
evenement etant produit en fonction, en outre, des 
intensites des photodetecteurs de deuxieme couronne. 

17. Procede selon 1 ' une des revendications 12 a 
16, l-'intensite du signal produit par chaque 
photodetecteur etant soit une integrale numerique du 
signal numerise, soit l'intensite au maximum du signal 
numerise, soit l'intensite au maximum du signal 
analogique . 

18. Dispositif pour le traitement de signaux 
produits par un ensemble de N photodetecteurs, 
comportant : 
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- des moyens pour comparer, pour chaque photodetecteur, 
l'intensite d ' un signal emis avec une valeur seuil, 

- des moyens pour comparer, pour chaque photodetecteur 
emettant un signal dont l'intensite depasse ledit 
seuil, l'intensite du signal qu'il emet avec 
l'intensite d'un signal emis par chacun des 
photodetecteurs voisins, 

- des moyens pour memoriser la valeur de l'intensite 
maximum parmi toutes les intensites, ainsi que les 
intensites des signaux emis par les photodetecteurs 
voisins du photodetecteur dont le signal presente 
l'intensite maximum-. 

19. Dispositif selon la revendication 18, 
comportant en outre des moyens pour produire un signal 
representat if d 1 une grandeur caracteristique d'un 
evenement a 1 1 aide des intensites . 

20. Dispositif selon la revendication 19, au 
moins un .signal- representatif de 1 * energie de 
1 ' evenement, ou de la position ou des coordonnees 
barycentriques de 1 1 evenement par rapport a 1' ensemble 
de N photodetecteurs, etant produit . 

21. Dispositif selon 1 1 une des revendications 
18 a 20, comportant en outre des moyens pour emettre un 
signal d' evenement, pour chaque photodetecteur, lorsque 
l'intensite est calculee. 

• 22. Dispositif selon la revendication 21, 
comportant des moyens pour detecter le passage du 
signal par son maximum et pour emettre le signal 
d ' evenement un certain temps apres le passage par le 
maximum. 

23. Dispositif selon 1 ' une des revendications 
18 a 22, comportant en outre des moyens de detection et 
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de re jet d'evenement saturants et/ou d 1 empilements 
d'ev^nements. 

24. Dispositif selon 1 ' une des revendications 

18 a 23, comportant des moyens de calcul et de 
memorisation des intensites des photodetecteurs de 
deuxieme couronne . 

25. Dispositif selon 1 1 une des revendications 

19 ou 20 et selon la revendication 24, le signal 
representatif d'une grandeur caracteristique de 
1 1 evenement etant produit a 1' aide, des intensites des 
photodetecteurs de premiere et deuxieme couronnes . 

26. Dispositif selon l'une des revendications 
18 a 24, comportant soit des moyens pour calculer une 
integrale numerique du signal emis par chaque 
photodetecteur, soit des moyens pour identifier 
l'intensite au maximum d'un signal numerique 
correspondant au signal emis par chaque photodetecteur, 
ou pour identifier l'intensite au maximum du signal 
analogique emis par chaque photodetecteur. 

27. Dispositif pour identifier une grandeur 
caracteristique d'un evenement se produisant devant un 
ensemble de photodetecteurs, comportant : 

— des moyens pour identifier un sous-ensemble de Ni 
photodetecteurs, ayant produit un signal en reponse a 
1 1 evenement, 

- des moyens pour produire un signal representatif de 
la grandeur caracteristique de 1' evenement, en 
fonction des intensites : 

* du signal d'intensite maximum parmi les Ni 
signaux, produits par les Ni photodetecteurs en 
reponse a l r evenement, 

* des signaux emis par les photodetecteurs voisins, 
dits de premiere couronne, du photodetecteur dont 
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le signal presente 1' intensite ■ parmi les Ni 
photodetecteurs . 

28. Dispositif selon la revendication 27, la 
grandeur caracteristique etant 1 1 energie de 
1 1 evenement . 

29. Dispositif selon la revendication 27, la 
grandeur caracteristique etant la position de 
1 1 evenement par rapport a 1' ensemble de N 
photodetecteurs. 

30. Dispositif selon la revendication '28, la 
position etant reperee par un couple de coordonriees 
barycentriques . 

31. Dispositif selon 1 1 une des revendications 
27 a 30, 1 1 intensite de chaque signal produit par 
chaque photodetecteur etant soit une integrale 
numerique de chaque signal, soit 1 1 intensite au maximum 
du signal numerique correspondant , soit 1' intensite au 
maximum du signal analogique. 

32. Dispositif selon l'une des revendications 
19 ou 27, les moyens pour produire un signal 
representatif d'une grandeur caracteristique etant 
propres a chaque photodetecteur. 

33. Dispositif selon l'une des revendications 
19 ou 27, les moyens pour produire un signal 
representatif d'une grandeur caracteristique etant 
communs a 1' ensemble des photodetecteurs. 

34. Dispositif pour la detection d'evenements 
comportant : 

— un reseau bi-dimensionnel de N photodetecteurs, 

- un dispositif selon l'une des revendications 18 a 33. 

35. Camera comportant un dispositif pour la 
detection d'evenements selon la revendication 34, les 
photodetecteurs etant des photomultiplicateur s . 
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36. Dispositif d'imagerie en rayonnement gamma, 
comportant une camera selon la revendication 35. 

37. Dispositif selon 1 1 une des revendications 
18 a 33, les fonctions de detection de 1 ! intensite du 
signal emis par chaque photodetecteur , de comparaison 
avec les photodetecteurs voisins, de memorisation et/ou 
de production d'un signal representatif d'un evenement 
etant assurees par un ASIC associe a chaque 
photodetecteur . 
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